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Abstract. We reconsider the electrostatic contribution to the persistence length, £e, of a single, infinitely 
long charged polymer in the presence of screening. A Gaussian variational method is employed, taking £e 
as the only variational parameter. For weakly charged and flexible chains, crumpling occurs at small length 
scales because conformational fluctuations overcome electrostatic repulsion. The electrostatic persistence 
length depends on the square of the screening length, £e ~ k~^, as first argued by Khokhlov and Khacha- 
turian by applying the Odijk-Skolnick-Fixman (OSF) theory to a string of crumpled blobs. We compare 
our approach to previous theoretical works (including variational formulations) and show that the result 
£e ~ K~^ found by several authors comes from the improper use of a cutoff at small length scales. For 
highly charged and stiff chains, crumpling does not occur ; here we recover the OSF result and validate 
the perturbative calculation for slightly bent rods. 



PACS. 36.20.-r Macromolecules and polymer molecules 
87.15.-v Molecular biophysics 



82.70.-y Disperse systems; complex fluids 



1 Introduction 

. Charged polymer chains, also called polyelectrolytes, have 
' many industrial applications (flocculants, viscosifiers, ad- 
I sorbants, etc.) related to their large solubility in water, 
I which is present even for a highly hydrophobic backbone. 
■ They are also important in biology since many biopoly- 
, mers (such as nucleic acids) are charged. 

In this paper we reconsider the controversial issue of 
] the dependence of the persistence length of polyelectrolytes 
on electrostatic interactions and, consequently, on the salt 
concentration in water, which has been of central interest 
during the past few years m^BHl- This is an important 
question with a large number of experimental implica- 
tions, since for instance, many properties of biological and 
synthetic polymers depend on their stiffness and thus on 
the salt concentration in the solution. It is therefore an 
important control parameter. Most of the theoretical at- 
tempts have treated electrostatic interactions on the linear 
Debye-Hiickel level, leaving aside the very important ef- 
fects connected with counterion condensation (which can 
be obtained on the mean-field level by properly taking 
non-linear effects into account) and with strong electro- 
static correlations, important when multivalent ions are 
present 0]|2||B]. Here we use the same simplifying assump- 
tions. Our purpose is therefore not to describe experiments 
or simulations (with explicit counterions) better than pre- 
vious calculations, but to detect and eliminate inherent 
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methodological shortcomings and inconsistencies in some 
of the previous calculations, which led to controversial re- 
sults. Our results should therefore not be compared with 
experiments, but rather with one of the many simula- 
tions of polymer chains where Debye-Hiickel repulsions are 
used jTIISllMTUnTT] . Especially the disagreement between a 
number of previous variational calculations is highly irri- 
tating |12II13II14II15II16II17II18| : variational calculations are 
typically known to be very robust, reliable and therefore 
useful in a wide range of physics applications. The imme- 
diate question therefore should be : is there any ingredient 
in fluctuating charged polymers that forbids the use of a 
variational theory ? The answer we give in this paper is 
clearly negative. Variational methods do work for charged 
polymers as they do for superconductivity, strongly cor- 
related electron systems, protein folding, random media, 
etc. This makes our methods in principle also applicable to 
more complicated problems such as polymer solutions and 
the coupling between counterion condensation and chain 
stiffening (which we will tackle in the future). 

Our calculational strategy consists in constructing the 
simplest variational Gaussian correlation function which 
i) reflects the stiffening of the chain at small length scales, 
but at the same time, ii) leads to a converging variational 
free energy and therefore makes introduction of a small- 
scale cutoff (as used in previous calculations) unnecessary. 
Quite surprisingly, the simple worm-like chain correlation 
function does not satisfy criterion ii) and therefore has 
to be discarded. The simplest workable correlation func- 
tion, related to the tangent-tangent correlation function, 
shows a two-scale behaviour, in agreement with previous 



Manoel Manghi, Roland R. Netz: Variational theory for a single polyelectrolyte chain revisited 



perturbative calculations 2 . It explains why the electro- 
static persistence length is difficult to extract in experi- 
ments or simulations. As our main result, we obtain that, 
in the flexible polymer case, where the chain crumples at 
small length scales and forms a persistent chain of Gaus- 
sian blobs (Gaussian-persistent regime), the electrostatic 
persistence length depends on the square of the screen- 



ing length, (e 



We thus confirm the scaling the- 



ory by Khokhlov and Khachaturian j^ who applied the 
Odijk-Skolnick-Fixman (OSF) ^(MT\ theory to a string 
of crumpled blobs. This is in accordance with previous 
theoretical studies [171111^ and very recent simulations [HI 
irnj . The uncontrolled usage of a cutoff at small length 
scales seems to be the origin of the very different results 
for the persistence length found in previous papers. In the 
next section, we summarize previous calculations for the 
persistence length of charged polymers. In Section 3, we 
present the general scheme of our calculations both for 
flexible and stiff chains. In Section 4, we study the case of 
flexible chains, using a variational Hamiltonian which ex- 
hibits a crossover from crumpled statistics at small scales 
to rod-like behaviour at intermediate scales. We show both 
numerically and analytically in the asymptotic limit of 
low ionic strength that ie ~ k~^. In Section 5, we re- 
turn to the case of stiff polymers, the persistent regime, 
and show why the worm-like-chain model cannot be used 
per-se for a variational approach. This is related to the 
behaviour at small scales. We use a slightly more compli- 
cated variational kernel which regularizes the small-scale 
behaviour and recover variationally the classical Odijk- 
Skolnick-Fixman result. Section 6 is devoted to the dis- 
cussion of our results and the comparison with previous 
theoretical calculations. Finally in Section 7, we give our 
concluding remarks. 



2 Previous Results 

The first studies of the influence of electrostatic interac- 
tions on polyelectrolyte stiffness have been done by Odijk pO] 
and Skolnick and Fixman 21. By performing a perturba- 
tion calculation on a slightly bent rigid charged rod using 
the Debye-Hiickel approximation, they found that the to- 
tal persistence length is 



i.e. for stiff polymers with a large bare persistence length, 
£b^o/^^ ^ Ij a-iid sufficient screening. In a systematic 
derivation, Barrat and Joanny found that the electrostatic 
stiffening is actually scale-dependent and ie as given by 
Eq. Q is only realized at large scales [T^. 

The case of flexible polymers, however, remained some- 
what unclear and a host of controversial results can be 
found in the literature. On one hand, Khokhlov and Khacha- 
turian J2| (KK) assumed that the OSF theory can be 
applied to a chain of polyelectrolyte blobs in the case of 
weakly charged flexible polyelectrolytes. At short length 
scales, the electrostatic repulsion is weaker than the chain 
entropy. The electrostatic blob size, ^ ~ Iqu^^'^ where n is 
the number of monomers (each of size £q) in the blob, is 
defined by the requirement that the electrostatic energy of 
the Gaussian polymer coil is equal to the thermal energy, 
that is {nio/AyiB/S, — 1. The blob size follows as 



A^ 



-1/3 



and the number of monomer in the blob is 

-2/3 



A2 



(3) 



(4) 



At a scaling level, the effective distance between charges 
along the direction of the string of blobs is then renormal- 
ized by the linear density of monomers, A'~^ ~ ^^A~^. 
By replacing A by A' in Eq. (01 , KK found 



Cp_ - 



KK 



A^ 



1/3 



io^^^ 



(5) 



ep = eo + 4 

where £q is the bare persistence length and 
static contribution: 



(1) 



the electro- 



t-OSF 



4^2 k2 



(2) 



The length is = e'^ / {AnekBT) is the Bjerrum length (the 
distance between two elementary charges interacting with 
the thermal energy ksT), A the distance between charges 
along the chain and k = [Att(.b{c+ + C-)]"'^/^ the Debye- 
Hiickel parameter, related to micro-ion concentrations, c± . 
This result is valid for polymer conformations which do 
not deviate too much from the rod-like reference state, 



which exhibits the same dependence in k^^ ^s the OSF 
result. This has been confirmed by Li and Witten who 
studied the effect of chain fluctuations in the presence of 
screening and showed that they do not affect the value of 
4 found by KK [ini. 

On the other hand, using a variational approach by 
modelling the flexible polyelectrolyte as a chain under 
tension, Barrat and Joanny |13j found a different n de- 
pendence: 

4 = ^BJ =^ - (6) 

K 

This result is consistent with a scaling argument proposed 
by de Gennes et al. [22] for solutions and has also been 
found by Pfeuty j23| using the Renormalization Group 
theory. More recent renormalization calculations have been 
done by Liverpool and Stapper [2^ and point to a simi- 
lar (even sublinear) dependence of the persistence length 
on K. A number of different approaches led to the same 
behaviour : i) an approach based on minimization of an 
approximated free energy jl2| . ii) a 1/d expansion (where 
d is the space dimension) j25) . and iii) variational ap- 
proaches 16 using the "uniform expansion" method first 
derived by Edwards and Singh J25II27| . 

It is important to mention that the result Eq. © poses 
problems at the crossover charge £_b^o — A^ , between the 
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Gaussian-persistent and the persistent regime. The k de- 
pendence for £e is then not continuous and a crossover 
formula should be found in theories. Among theoretical 
works which predict Eq. JHJ, there is no satisfying answer 
to this question. 

Recently, Netz and Orland J7] used the most general 
Gaussian kernel for the variational Hamiltonian |2H| and 
recovered the well-known three distinct scaling regimes, 
defined according to the value of the screening parame- 
ter k£o and the charge parameter iB^o/A^- For weakly 
charged polymers and at large screening, one is in the 
Gaussian regime where the notion of electrostatic persis- 
tence length is meaningless. At small screening, and ac- 
cording to the charge and the stiffness of the chain (for 
(,b(-o/-^ > 1); one reaches the persistent regime where (.e 
is given by the OSF result, equation (0). The Gaussian- 
persistent regime is reached when the electrostatic repul- 
sion between monomers is weaker (for iB^o/A^ < 1) and 
the chain is crumpled at small length scales. It is shown 
that in this last regime, the asymptotic electrostatic per- 
sistence length (when k(.q — *■ 0) is given by the KK result, 
equation ©. 

Although the work by Netz and Orland is in accord 
with the growing consensus that the KK result is asymp- 
totically correct, it uses quite complicated mathematics 
and cannot be easily extended to other systems, such as 
polyelectrolyte solutions. Moreover, the persistence length 
is found quite indirectly by matching the small-distance 
behaviour with the asymptotic swelling range due to ex- 
cluded volume interactions. Hence, the intermediate dis- 
tance range where the chain statistics is truly Gaussian 
at scales larger than the persistence length is not consid- 
ered. However, such an intermediate range always exists 
provided that £p 3> k~^ which is typically the case both 
in the persistent and the Gaussian-persistent regimes [3]. 
The persistence length should thus be determined by the 
crossover between the persistent range and the interme- 
diate Gaussian range. Finally, in order to compare to ex- 
periments and simulations, it would be interesting to find 
the K-dependence also near the crossover, k^ ~ 1, where 
screening becomes important, and not only asymptoti- 
cally. These are the reasons that motivated this study 
where the variations of ^e(K) are determined self-consistently 
by choosing (.e as the variational parameter. 

In a number of papers a similar single-parameter vari- 
ational theory was used, but the variational persistence 
length refiected the swelling behaviour of the chain for 
large scales 29,30,31,32,33,34 . This is borne out by the 
fact that the persistence length depends on the chain length 
with a characteristic power law £e ~ N"* which in the 
infinite chain limit TV — > oo gives the correct swelling 
behaviour for the radius of gyration. In these calcula- 
tions, the persistence length does not reflect the mechani- 
cal properties of the chain at small length scales, because 
one would expect the persistence length to saturate to a 
finite value as TV — > oo. In our paper, in contrast, we fix 
the large-scale behaviour as being Gaussian (thus corre- 
sponding to the actual intermediate regime for real chains) 
and the variational parameter reflects the transition from 



the rod-like to the Gaussian behaviour. The reason for 
this choice is that taking into account both the stiffen- 
ing (at small length scales) and the swelling (at large 
length scales) requires two distinct variational parameters, 
a formidable task which can be done only asymptotically 
and within uncontrolled approximations |17|. 



3 Model 

In the following, we consider a Gaussian charged poly- 
mer in a solution of monovalent salt (concentrations c-|_ 
for cations and c_ for anions) . The polymer chain of poly- 
merization index TV and a chain length L = N£q is sup- 
posed uniformly charged with quenched charges +e. The 
distance between charges along the chain is given by A. 
We remain at the linear level for electrostatic interactions, 
assuming a Debye-Hiickel screened interaction between 
charges due to the presence of monovalent counterions 
and salt. As explained in the Introduction, we therefore 
neglect non-linear effects connected to counterion conden- 
sation at low salt concentrations (we plan to incorporate 
such effects in the future by using a modified variational 
action; in the present paper we mostly try to resolve in- 
consistencies in previous treatment of this linear model) . 
We use a continuous description of the polymer, param- 
eterized by the curvilinear index s and where r(s) is the 
position of the monomer labelled by s. The Hamiltonian 
of the system is ,2^-2^211 



pn[r] 



PH,i[r] + f3Hstis[r 

1 '•^ - 

2 



ds 



/o 
2A2 



/5Hi„t[r] 



i-/M + '^-^'W 



Jo 



\v{s) 



r s 



where [3 = {kBT)~^. The first term is the standard elastic- 
ity of entropic origin, the second one the bending stiffness 
of the chain, and the last term the electrostatic interac- 
tion. The partition function is 



Z = 



fvrexp{-(3H[r]} 



(8) 



and is intractable due to the electrostatic term. 

To make progress, we choose a variational Gaussian 
Hamiltonian which most generally reads 

f3no[r] = l f ds f ds'r{s)g-\s^sy{s') (9) 
^ Jo Jo 

where g~^ is the Gaussian kernel. Gaussian statistics im- 
plies that g is the monomer-monomer correlation function: 



(r(s)r(s'))o = 3.9(,s-s') 



(10) 



where (...)o denotes the expectation value computed with 
the variational Hamiltonian J^. It is related to the squared 
monomer separation as 

{[r{s) - r{s')]\ = G{s - s') = 6.g(0) - 6g{s - s') (11) 
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The tangent-tangent correlation function is simply 



(r(s)f (.'))o = 3 



32 



1 



dsd. 



l9{s~s') = -G"{s-s') 



(12) 



weaker than the thermal energy, /csT, for {t^ / A)^ £ b / ^o < 
1, and thus is not sufficient to align the two segments. 
This threshold is equivalent to the condition 



The variational free energy reads in the Gibbs-Bogoliubc 
form 



^2 



< 1 



(17) 



F„ 



Fo + {H-Ho)o 



(13) 



where (3Fo = — In Zq is the free energy associated with 
the variational Hamiltonian defined in equation @. The 
Gibbs inequality ensures that Fvar > -Foxact when i^var is 
minimized with respect to the variational parameters. As 
we are interested in the persistence length, we will choose 
the electrostatic contribution £e (to be defined below) as 
the only variational parameter. The value of the persis- 
tence length will thus be determined by minimizing equa- 
tion (|13|l with respect to £e- It must be emphasized, how- 
ever, that we will restrict our choice of variational Hamil- 
tonians to the subclass of non-swollen (ideal) correlations 
at large scales, and by doing so, we do not obtain the true 
minimum of the variational free energy. 

Before starting explicit calculations, it is useful to con- 
sider the Hamiltonian Q and to rescale all lengths by the 
bare persistence length according to r = £of and s = £os. 
We obtain 



/3H[r] = 



3 .N 

2 A^ 



ds 



i2 -^2 ■ 

r (s) +r (s) 



(14) 



N 



dS 



N 



, exp(-K4|f(g)-f(gO|) 
|r(S)-r(S')| 



with three dimensionless parameters, namely the charge 
parameter, £b£o/A'^, the screening parameter k£o and the 
polymerization index N. In this paper, we consider the 
thermodynamic limit N —^ 00 and therefore are only left 
with two parameters. 

In the limit £3(^0/^'^ — the tangent-tangent corre- 
lation function can be calculated exactly from Hamilto- 
nian Q and reads 



(15) 



(r(s)f(0))wic = exp ( - — 



while the squared monomer separation is given by 



Gwic(s) = 2£l 



^-1. 



""P'~4 



(16) 



The correlation function therefore is identical to the worm- 
like-chain (WLC) model which has been first proposed 
by Kratky and Porod for describing neutral semiflexible 
chains "3^1 . 

In our variational calculation, we will choose a mod- 
ified WLC Hamiltonian, or, equivalently, a modified cor- 
relation function G(s), which will reflect non-zero values 
of the charge parameter £b£q/A'^. The choice of the vari- 
ational correlation function will contain our knowledge of 
the behaviour of charged chains. For weakly charged and 
flexible chains, the electrostatic repulsion between two ad- 
jacent charged segments of length £0 and charge £q/A is 



In this case, the statistics of the chain is Gaussian at small 
scales and the chain has locally a crumpled configuration. 
The electrostatic blob size |211, £,, is thus defined at the 
scaling level by Eq. Q. This blob size is found using a 
scaling argument but can also be found variationally as 
shown in Refs. |22lll3j plus a logarithmic correction jl7| . 
In the following, we neglect logarithmic corrections and 
use Eq. (@J) for n in the rest of the paper. For weakly 
charged chains, we choose a variational correlation func- 
tion in accordance with these scaling results and which 
exhibits the crumpled configuration at small scales (Sec- 
tion 4). The case of highly charged and stiff chains for 
£b£o/A'^ > 1, which is treated in Section 5, is more subtle 
since here the variational correlation function has to be 
treated in a way such that all integrals occurring in the 
variational free energy are regular in the small distance 
limit. 



4 Flexible polyelectrolytes 

At scales smaller than the electrostatic blob size (s < n) 
the chain conformation is Gaussian, as argued above. For 
larger scales {s > n), the chain shows, for sufficiently large 
charge parameter, electrostatic stiffness and should be de- 
scribed by an analogue of Eqs. ifT^ and lfTC|l . provided 
that we properly rescale the contour length by the blob 
size. 

To simplify the mathematics, we in this section use a 
flexible polymer model. 



^^W^if^^^^ 



(s) 



-^ f ds 
2^ Jo Jo 



'^^., exp(-^|r(g)-r(g')|) ^^g^ 



|r(s)-r(s')| 



The simplest variational tangent-tangent correlation func- 
tion can be written as 



1 



(f(s)f(0))o = -G"(s) = 4-5(5)- 



1 



■ exp 



l-s| 



1I/24 



(19) 



which reproduces the crumpled statistics following from 
Eq. (|18|l at small scales, and also a persistent behaviour 
with a persistence length £e (or n^^'^£e along the contour 
of the chain) at large scales. The factor n~^ in front of 
the second term ensures the correct crossover between the 
Gaussian and rod-like behaviour. By integrating twice, we 
obtain the mean-squared monomer-monomer separation 



G{s)=£o\s\ + 2£l 



i^/H, 



— 1 + exp 



11/2^ 



(20) 
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We remark that our variational choice is exactly the same 
as the one proposed by Barrat and Joanny IS . This can 
be shown explicitly by integrating out the fluctuating ten- 
sion t(s) (which can be done exactly) in their partition 
function. 

In the following, we compute the different terms of the 
variational free energy equation IjlSI) using equations 10 , lU 
and p9|l . We denote by h{uj) the Fourier transform of the 
tangent-tangent correlation function 



By inserting equation H19|) in the last equation, using the 
definition of n and after a few algebraic calculations we 
obtain 



N 



/3(Hi, 




h{Lo) = I ds(f(s)r(0))oe'"' 
dsG"{s) cos(a;s) 



where 



where we used that G" is an even function. We find 

e, 2 



Hlo) = io 



(21) 



(22) 



C 1 + Ili^UJ^ 

The average value of the elastic part is easily found to be 
/3(Hc:)o= rf^|^M-) = |^(r(0)r(0))o 



3A^ 

~2~ 



duj 3N 
27r 2~n 



(23) 



which does not depend on the variational parameter. Fur- 
thermore, it is easy to check that (7io)o is a constant in- 
dependent of ie- The entropy term is given by 



2 ./ _ 27r 



h{u;) 



niU^ 



(24) 



By substracting the free energy part associated with the 
ultraviolet divergence (w -^ oo), i.e. the small length-scale 
fluctuations which do not depend on the variational pa- 
rameter, 

/3F,. = ^ r ^Inielco^) (25) 



2 7-00 27r 



we find 

P [Fo - F^ 



3N ^ 
2~nTe 



27r V O + a;2 




(26) 



In the limit N ^^ oo, the interaction term reduces to 



P{'Hint)o 



LIb 

Ub 
A^ 



As 



As 



exp(— K|r(s) — r(O)l) 
|r(s)-r(0)| 
d^q 47r 



(27) 







(27r)3 K^ -h q- 



2 _i_ri2 

2 



^jq|r(s)-r(0) 



«b2 p, Z-"", q 



e-5«'G(s) 



G{x) ^ xf + 2{e-^ ~ 1 + x) 



and the complementary error function is 

2 f°° 
erfc(x) = —;= / diexp(— i^) 

y/T^ Jx 



(29) 



(30) 



At scales larger than the electrostatic blob size ^, tangents 
are correlated over distances up to the persistence length 
£e if the Debye length, k~^ , is larger than ^, that is, if 
K^ < 1. The case where k^ > 1 and n > 1 corresponds to 
the Gaussian regime ^21 and the notion of electrostatic 
persistence length is then meaningless. We have checked 
that, indeed, t^ goes to for k^ > 1. 

The minimization of the variational free energy per 
blob, nPF^i^r/N, with respect to £e is first done numeri- 
cally. The result £e{i^) for k^ varying from 0.001 to 1 is 
shown in Figure Q in a log- log plot. We observe that the 
electrostatic persistence length decreases with k following 



the power law If, 



The solid line denotes the function 



in the limit of small k^ (k^ < 0.1) 



(31) 



which demonstrates that the result found by Khokhlov 
and Khachaturian by scaling arguments is correct for De- 
bye lengths larger than 10 blob sizes. For smaller De- 
bye lengths (0.1 < k^ < 1), our results show a smooth 
crossover to the results obtained in Ref. TT, Eq. ^, 
shown in Figure^as a broken line. It should be noted that 
this part of the curve (0.1 < k^ < 1) varies depending on 
the variational choice and is not universal (see Appendix). 
To understand where the discrepancy between our asymp- 
totic scaling and the result by Barrat and Joanny comes 
from, we analyze our variational free energy in the limit 
of weak screening. 

From Eq. H26(l . we can approximate the entropy per 
electrostatic blob by 



^«F. 



i^uv) 



V2 Ve 



1/2 



(32) 



in the limit ie/^ ~* oo, which is the low-screening limit. 
Obviously, this limiting formula is different from the re- 
sult found by Barrat and Joanny for the entropy, ^(3Fq ~ 
\n{ij^) [see Eq. (7) of Ref. [fS]]. This different depen- 
dence of the entropy on l^ is the sole reason for their differ- 
ent scaling of the persistence length, resulting in l^ ~ k~^. 
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(36) 



where the structure factor, Sg^{q), of an infinite chain 
computed using the variational Hamiltonian, TYq, is given 

by 

' jq|r(s)-r(0) 

' 



5o-(g) = 2^"ds(e 



(37) 



K^ 



Fig. 1. Numerical minimization of the variational equa- 
tion 1131 1 with respect to le'. log-log plot of the renormalized 
persistence length, £e/^, vs. the renormalized Debye parame- 
ter, K^ (where ^ is the electrostatic blob size). The solid line 
denotes the function £e = 7^ while the broken line is propor- 
tional to k,~^ . 



In the limit £e/S, —^ 00, the chain can be viewed as a semi- 
flexible chain made up of N/n rods of length ^ and of 
persistence length i^- We thus neglect the internal blob 
structure. We can use an interpolating formula for the 
polyelectrolyte structure factor, S^{q) = (1 + qie)/{q^ie) 
which exhibits the correct limiting behaviour for small mo- 
menta in l/{q'^£e) and for large momenta in 1/g PSj. Using 
a upper cutoff at g ~ ^^^ in Eq. (|36|) . we easily find 



N 



m 



1 



int/O 



HO 



TT 1 
4^ 



o 



1 



(k4 



(38) 



Taking a derivative with respect to ^e and inserting the 
KK scaling Ea. (|31|l . we recover Eq. H34|l . Hence, we note 



that the 



The interaction term equation H28|) is more difficult to 
evaluate. The two terms cannot be computed separately 
because of a cancellation of divergences. To make analyti- 
cal progress, we assume the result (|5T|l to be valid and re- 
place ^e by 1/{^K?) in the first two terms of equation (|28l) . 
We then look at the limit k^ — > 0. We find numerically for 
the leading term 



dependence of te comes directly from the 



balance between the entropy of the chain in equation (|32|1 
and the electrostatic energy in equation l|38|l in the limit 
of weak screening. 



^^mHjidi, 



-^/3(Wint)o 



HO 



(33) 



It corresponds to the electrostatic energy of the blob in a 
cylindrical geometry and diverges logarithmically when k 
goes to 0. This term is independent of the electrostatic per- 
sistent length. To find the next leading term, we calculate 
the asymptotic behavior of the derivative of -^/3(7iint(^e, k))o 
with respect to £e and then assume the KK law. As is vi- 
sualized in Figure 13 we find 




«_^^ ^(^mt)o 



-«)^ + 0(«)5) (34) 



--{K^ir 



Keeping the first terms of the asymptotic expansion of 
dFvai/d£e, we thus get 



dFh 



d{H 



int/O 



dt. 



^1/2 
772 



k3^2 ^ (35) 

and thus recover Eq. H31|l in the asymptotic limit k^ -^ 
(within a numerical coefficient of the order of unity) in a 
self-consistent way. 

Besides this asymptotic but cumbersome way, it is also 
possible to find the leading term in (/t^e)"^ for the inter- 
action term by using an approximated structure factor for 
the semi- flexible chain ^31 ■ Starting from Eq. H27|l , we see 
that 



/3(Hi, 



Ub 



ds 



d^q 47r 

(27r)3 k2 + c 



,jq|r(s)-r(0) 



Fig. 2. Numerical check of Eq. 1341 . The first deriva- 
tive of the interaction term per electrostatic blob, 
■^/?5c'(7iint)o/9^e|f^=(K25)-i, decreases as —{i^iY when 
le is assumed to follow the KK law, (.^ — ('^^C)~^- 



5 Stiff polyelectrolytes 

In the persistent regime, for stiff and strongly charged 

polymers, £o-^bM^ > 1, and at small screening, £o-^bM^ > 
(^0/^)^5 the above variational kernel must be modified. Ac- 
cording to OSF theory, which is based on the electrostatic 
contribution to the energy of a uniformly bent rod, the 
electrostatic monomer-monomer interactions are already 
relevant on length scales comparable to the bare persis- 
tent length, leading to an effective persistent length given 
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by equation ^. However, Barrat and Joanny have shown 
within a quadratic perturbation treatment that the persis- 
tence length is in fact scale-dependent and that the OSF 
prediction corresponds to the large-scale limit ^Hl- Since 
all these calculations are perturbative and only valid for 
weakly bent chains, it is interesting to have a variational 
determination of the effective persistent length £p = £o+ie 
in this persistent regime. The first choice would be the 
WLC model with a variable persistence length. Treat- 
ing the persistence length £p as a variational parameter, 
the tangent-tangent correlation function and the squared 
monomer separation are given by Eqs. H15(l and H16(l where 
£o is replaced by ip. For an infinite chain, A^ ^ oo, a Gaus- 
sian Hamiltonian which leads to these equations is Eq. |(7J 
with Ib/A-'^ = and £p instead of Iq. However, this vari- 
ational Hamiltonian with £p as the variational parameter 
does not work because the variational entropic free energy 
diverges at small length scales 



/3-Fwic — — 



3L 



duj 

2^ 



In 



^^^^2(1 



■^2^2 



(39) 



since the integrand behaves like ln(^o^p<^'*) when w -^ oo. 
This divergence depends on the variational parameter, and 
therefore, a variational approach based on such a Hamil- 
tonian is ill-defined. This fact is reflected by Monte-Carlo 
simulations and perturbative calculations for the tangent- 
tangent correlation function of a charged chain : at small 
length scales, the decay of this correlation function is given 
by the bare (mechanical) persistence length, and only at 
larger length scales one finds a crossover to the electro- 
static persistence ^3]. This can be described by a scale- 
dependent persistence length £p{uj) = £o + £e{uj) where 
£ei^) — for oj — > oo and ^e(^) — ^osf for uj —^ which 
leads to a tangent-tangent correlation with the limiting 
behaviours 



1 



s| 



for \s\ < Sc 

(f (.)f (0)) ^ <; f |,| ^ ' ' (40) 

where the crossover contour length s^ has been determined 
in Ref. 13 and is 



1 



k^/T+£osf/£o 



(41) 



It transpires that the variational correlation function has 
to be chosen in accord with these limits. A simple and con- 
tinuous choice for the tangent-tangent correlation function 
which fulfills Eq. ^ is 



(f(s)f(0))o = Bexp 



4 + 4 



(42) 



(l-B)exp 



£e + {l-B)£o 

hoiie + £a){l - B) 



where £e is the variational parameter and the factor B is 
chosen such as to satisfy continuity 



B = exp 



«;V4(4 + 4)3/2 



(43) 



Hence, the crossover Sc — V^/i'^V^e + 4) depends on 
the variational parameter. The case of a neutral semi- 
flexible polymer is recovered by taking Sc —^ oo (which 
implies S = 0) in Eq. 1)42(1 . The tangent-tangent corre- 
lation function is plotted in Figure 13 for k£o = 1 and 

4 = 104- 



<i-(5)r(0)) 
1 




Fig. 3. Log plot of the variational tangent-tangent correlation 
function vs. s = s/£q [see Eq. H42^ ]. The broken lines are the 
two asymptotic behaviours given by Eq. MUll . The parameter 
values are hIq = 1 and ie = 10 



The Fourier transform, defined in Eq. (121(1 . is easily 
computed 

4 + 4 



h{u;) = 2B 



1 + (4 + 4)2^2 

2(1-5)2(4+4) 



(44) 



4 



4 + (1 - B)4 



1 



(£e+eo){i-B)eo 
e„+{i-B)eo 



Following the lines of Section 3, we obtain for the entropy 



l3iFo - ^„) 
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i + 4w2)/i(w) 



37V 



2(4 + 1) 



1-B 



(45) 



(i-B + 4)(i-S + B4 



V4+T(i - B) 

where £e is the adimensional electrostatic persistence length 

4 = ^ 
4 

and 



(46) 



PFn 



3L 



duj 
2^ 



In 



1 



u;'£f,{l+£f,Lj^) 



(47) 



is the entropy of a neutral semi-flexible polymer. The elas- 
tic term is 

/3(Wei)o = ^(^-4^)(f(0)f-(0))o 



3iV 



i?£? 



4 (1-B)(4 + 1)2 



(48) 
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The electrostatic energy is found using Eq. 1)27(1 which used in the simulations by Barrat and Joanny. An expla- 
yields nation could be that close to the crossover Ib^q/A^ = 1, 

the correction, cp, to the OSF law according to 

p p poc I /"o" 1 



X erfc 



k(4 + 1) 
V3 



vm 



4 A2 (k4)2 



A^ 



\K'(£,+iyfix) 



<49) 



where 



fix) =3(0-=" -1 + x) 



(l-Bf 



(4 + 1-5)2 



1-B 



1 -B 



- 1 



The interaction term, Eq. (|49l) . diverges logarithmically 
for small a; (rod- like divergence). In the computation, we 



rCK 

Jo 



dx 



thus subtract N^y- j^ j^^^p^ 
pend on the variational parameter. 



looa 



which does not de- 



^/L 




Fig. 4. Numerical minimization of the variational equation 
free energy with respect to ie for stiff chains : log-log plot of 
the renormalized persistence length, £e/£o, vs. the renormalized 
Debye parameter, kIq for two values of Ifjio/A^ : 2.25 (open 
circles) and 4 (full circles). Solid lines are the power law k^^. 



The minimization of the variational free energy is done 
numerically and the result is plotted in Figure^for €b^o/^^ 
2.25 and 4. We observe that the OSF result is recovered 
with a decrease i^ ~ k~^ in both cases. 

We should add that we find the persistence length to 

depend on the charge parameter €^4/^^ asfg ~ (^^4/^^) 
which is different from the linear dependence found by 
OSF [Eq. ©I- At this point we have no clear interpre- 
tation for this discrepancy. Because of numerical limita- 
tions, we could calculate the effective persistence length 
only for moderate values of the charge parameter up to 
£b£oM^ = 8, which is the range of parameters found for 
fully charged flexible (synthetic) polymers. This parame- 
ter range is thus very far from the value £3^0/ A'^ = 150 



(51) 

is very large. This point will be studied by suitable simu- 
lations in the future. 



(50) 5 Discussion 



In this section, we discuss the results of Section 3 and 
compare our approach to those of Barrat and Joanny |13| , 
Schmidt ^^, Ha and Thirumalai 16,18 and Muthuku- 
mar et al. |29ll3()j . All these works are variational ap- 
proaches with only one variational parameter. However, 
as explained in the Introduction, we shall distinguish be- 
tween two types of works, according to the physical inter- 
pretation of this parameter. In Ref. |29II30| . it reflects the 
swelling behaviour at large length scales and reproduce 
the correct swelling law as a function of A^. This varia- 
tional parameter is then not interpreted as the true me- 
chanical persistence length. In the other references |12II13I 
I1()II18| . the variational parameter is interpreted as a per- 
sistence length and the swelling behaviour at large scales 
is not considered. Our variational calculation belongs to 
this category. We show that, in this case, the persistence 
length is finite for an infinite chain, A'^ — > 00. 

As mentioned in Section 3, our variational calculation 
is similar to Barrat and Joanny calculation. They find the 
result Eq. 1(38(1 for the interaction term, by approximat- 
ing the structure factor of a polyelectrolyte and by using 
an upper cutoff at g ~ ^~^. Hence they "smooth" elec- 
trostatic interactions inside the electrostatic blobs. This 
approximation is valid only in the limit 4/C ~* OO; i-C- 
in the limit of weak screening. However, using the same 
cutoff at (7 ~ ^^^ in the calculation of the entropic term, 
they find /3Fo ~ — ln(£e/^), which is different from our re- 



sult, Eq. l|^. and leads to their law 



^ Evidently, 



this un-controUed cutoff leads to the correct result for the 
interaction term but not for the entropy. 

Schmidt (12j use a different variational formulation, 
where the elastic energy is calculated for a stretched ran- 
dom flight chain and the electrostatic energy is approx- 
imated using an hybrid rod-coil model. This leads to a 
weaker dependence than k~2 foj- £^ but which does not 
follow a power law and tends to vary as k~^ at high ionic 
strength. However, l^ increases with the contour length 
^ and seems to diverge for TV —> 00. Moreover, as carefully 
explained in Ref. 12 , this calculation does not apply to 
the low ionic strength limit in view of the approximations 
made. 

Ha and Thirumalai [16II18( use the same variational 
squared monomer separation, Eq. ((20() . but follow the uni- 
form expansion method, first developed by Edwards and 
Singh (26II27( . This perturbative variational method is per- 
formed on the mean-square end-to-end distance, (i?^). It 
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consists in calculating self-consistently the persistence length 
by requiring (i?^) — (i?^)o to be zero at first order in 
^o^s/^^- This method has proved to be powerful in the 
calculation of the Flory exponent for the end-to-end dis- 
tance of neutral polymers. For our issue, Yethiraj shows 
that in the limit L/l^. -^ oo the Gibbs-Bogoliubov ap- 
proach and the uniform expansion method applied to the 
end-to-end distance lead to the same self-consistent equa- 
tion |31| . Clearly, both variational approaches lead to sim- 
ilar results for a given observable. The important point is 
to carefully consider, in both approaches, the behaviour 
at small length scales. 

This method is applied both in the limit £o^b/^^ ^ 1 
and at the so-called non-asymptotic limit, fb^oM^ ~ 1, 
when the blob size is of the order of the bare persistent 
length ^ ~ €o 16,18 . In the first case, they find the KK 
result (18) but in their calculation they approximate the 
chain by a rod on a scale of the order of £e (thus factor- 
izing the statistical weights, see Eq. (2.18) in Ref. 1181). 
Moreover they skip a priori the integrals on long length 
scales, s > i'e, in the calculation, arguing that they lead 
to an effective excluded- volume interaction. Hence, the de- 
pendence of the persistence length on N is eliminated. In 
the non-asymptotic limit. Ha and Thirumalai use a simple 
Gaussian variational weight used for neutral and flexible 
polymers and find £e — (^o-^bM^)^^^'*^^- However, as we 
explained in the preceding section, this choice of Hamil- 
tonian leads to a divergent entropic term. To circumvent 
this problem, the evaluation of the statistical averages is 
done by imposing a lower cutoff at q ~ i~^ . We believe 
that this approximation leads to the scaling i^ ~ k~^. 

As shown by Muthukumar et al. |29II3U) the uniform 
expansion method applied to (i?^) cannot be used to infer 
the scaling of the persistent length. Indeed, to make the 
variational equation tractable (see Eq. (3.38) of Ref. 29 ), 
they approximate the variational kernel h[ijj) (noted li{q) 
in their paper) by the effective persistence length i,. [valid 
in the limit cj — > 0, as it can be seen from Eq. H22|l ] and 
they evaluate the integrand for w ~ I-k jL. Hence, these 
approximations are valid in the limit of large length scales 
and this method is powerful to get the large scale be- 
haviour : they indeed recover the swelling exponent in the 
infinite limit N ^ oo, (i?^) ~ N^^^ due to effective ex- 
cluded volume interactions (which are the screened elec- 
trostatic interactions if k~^ is non zero). But the extrap- 
olated persistence length from (i?^) ~ ipN is ip ^ N^^^, 
and does not correspond to the mechanical persistence 
length. In this approach, the small scale behaviour is, in 
a sense, averaged out, and all the electrostatics contribute 
to excluded volume interactions. 

The uniform expansion can be used to calculate next- 
leading corrections in an asymptotic expansion of the per- 
sistence length ie of the order of 0((k^)~^. But in this 
case, it must be applied to to the observable conjugated 
to the variational parameter, ig, namely the bending en- 
ergy /dsr^(s), and not to the end-to-end distance. This 
could be valuable near the crossover, k^ ~ 1, between the 
Gaussian-persistent regime and the Gaussian regime. 



In two recent papers |??lll()j . simulations on long poly- 
electrolytes (up to TV = 4096) show a behaviour consis- 
tent with the KK formula (jsj. Discrimination between 
£e ~ K~^ and £e ^ k~^ is clearly seen only for very long 
polymers which is one of the reason why the available 
experimental data are not very conclusive on this point. 
Furthermore, in experiments there are other effects such 
as counterion condensation and electrostatic correlations 
which are not included on the Dcbyc-Hiickel level. 



7 Concluding remarks 

Our calculation is done within three important assump- 
tions : i) we study an infinitely long chain {N —^ oo), ii) 
we do not consider excluded volume interactions at large 
length scales and, iii) we take into account electrostatic 
interactions at the linear level [a Debye-Hiickel potential 
is used in equation 0]. Within these assumptions, we 
show that the electrostatic contribution to the persistence 
length, £e, is proportional to k~^, both for flexible (KK 
result), Eq. ©, and stiff polymers (OSF result), Eq. Q. 
This result is very robust (we tried a different variational 
choice shown in the Appendix which leads to the same for- 
mula) and is in contradiction with several theoretical pa- 
pers where a k~^ dependence was found for flexible poly- 
mers. We show that this difference can be explained by 
the use of a cutoff at small length scales in those works. 

Obviously, a full comparison with experimental mea- 
surements on very long chains would be necessary to val- 
idate this variational approach. Comparison with experi- 
mentally determined persistence lengths of polyelectrolytes 
has been done in several articles JEIEIEI- Essentially, 
some experimental papers raised some doubts about the 
validity of the OSF result and proposed a law in £e ~ 
K~^ [37] . This law has been found using approximate for- 
mulas for the radius of gyration to deduce the persistence 
length, and the discrepancy can be due to an incomplete 
theory used in the data analysis, as shown by Ghosh et 

al. [sni. 

More importantly, a deviation from the Debye-Hiickel- 
type theory used in our paper arises due to non-linear 
electrostatic effects. For highly charged polymers (espe- 
cially in the persistent regime), the renormalization of the 
bare charge of the chain associated with the condensation 
of counterions (known as the Manning condensation |38p 
also influences the dependence of £e on the Debye screen- 
ing length. A recent study of this effect using a varia- 
tional approach applied to a charged cylinder [^ shows 
that the effective charge density along the polymer fol- 
lows A ^ K^'^'^ which would lead to a modified power 
law £e ^ K~^'*. We are currently studying this effect for a 
semi-fiexible polyelectrolyte. 

Finally, we hope that this simple variational approach 
can be extended to the issue of the electrostatic persistent 
length of polyelectrolytes in semi-dilute solutions by com- 
bining this type of variational approach with the random- 
phase approximation to account for monomer-monomer 
correlations and screening induced by neighboring poly- 
mers 
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Appendix 

In this appendix, we show that the KK result is also found 
with a piece-wise continuous tangent-tangent correlation 
function which exhibits the same behaviour as the mod- 
ified WLC model at the asymptotic limits \s\ <^ n^/^^e 
and \s\ 3> n^/^^e- The advantage of this formulation is 
that the interaction term, Eq. (|28|l is easily computed, 
and therefore we can show that the large-scale term in 
the integrals is negligible in the determination of the per- 
sistence length. Moreover, we thus prove that our result 
for the salt-dependence of the persistence length is robust 
and insensitive to the detailed choice of the variational 
correlation function. 

Our variational choice for the squared monomer sepa- 
ration is 



(jpw{sj 




1 \s\ 
3 n^/^e^ 



for 
for 



|s| < ni/24 
|s| > ni/24 



where "pw" 



(52) 
refers to piece- wise. For \s\ < n^'^ie, the 
first line of Eq. (|52|l is a sum of a linear term (Gaussian 
behaviour at small length scales) and a term which im- 
poses almost a rod-like behaviour for arc lengths ranging 
in n£o < |s| < n^'^^e- These linear and quadratic terms 
in the first line of Eq. (|52|l are of equal magnitude for 
s = n£o, the number of monomers in a Gaussian blob at 
the onset of the persistent behaviour. For \s\ > n^^'^le, the 
statistics is again Gaussian, as expected for a semi-flexible 
chain at large length scales. 

The tangent-tangent correlation function is thus com- 
puted through (r(s)f (0))pw = ^Gpwi^) 



(f(s)f(0))f 



^m + ^ ( 



1 



^1 



n^/^l. 



-] for \s\ < ni/24 

for |s| > ni/24 
(53) 
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Fig. 5. Same plot as in Figure^with a piece-wise variational 
squared monomer separation function Gp„ 



The function h, defined in Eq. (|21(l . is then 
2 l-cos(tjni/24)' 



hpwiuj) = ia 



1 



4C 



(54) 



We find exactly the same elastic part, Eq. (|23|l . as in the 
modified WLC model and for the entropy term we have 



P i^-^pw ^uv) 



47r£p n 



dx In 1 



^e 1 ~ COS X 



where F^^ is given in Eq. (|25|l . The interaction term is 
found to be 



-^/3(Wint)pw 




The numerical minimization of Eq. H13|l is shown in Fig.|Sl 
and shows that the KK result is also found with this vari- 
ational approach. 
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Fig. 6. First derivative of the interaction term per electrostatic 
blob,-^/3^c'(7-^int)o/i9^e|f^=(re25)-i, vs. K^ when (.^ is set to fol- 
low the KK law, le — {i^^Cj~^ i^i ^ log-log scale. The dotted 
line is the power law (k^^- 



We have checked numerically that the entropy per elec- 
trostatic blob given by Eq. (|55|l has the same dependence 
in —\Ji,lie for £e/^ — > 00 as for the modified WLC model 
[Eq. IjSHjl ]. The last term of equation (|^ corresponding to 
the large scale behaviour, i.e. for \s\ > n^/^£e, can easily 
be evaluated using the saddle point expansion and varies 
like ~ (ftfe)"^- By applying exactly the same routine as 
in Section 3, we show (Fig. ^ 



^/3(Wint)0 ■ 



ln«).- 



(57) 
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which yields the same result as in Section 3. We thus note 
that the electrostatic energy at small length scales, \s\ < 
n^''^£e, is larger in the limit ie/^ -^ oo and governs the 
dependence of ig with k. 
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